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In two-dimensional systems with a continuous symmetry the Mermin-Wagner-Hohenberg theorem 
precludes spontaneous symmetry breaking and condensation at finite temperature. The Berezinskii- 
Kosterlitz-Thouless critical temperature marks the transition from a superfluid phase characterized 
by quasi-condensation and algebraic long-range order to a normal phase, where vortex proliferation 
completely destroys superfluidity. As opposed to conventional off-diagonal long-range order typical 
of three-dimensional superfluid systems, algebraic long-range order is driven by quantum and ther¬ 
mal fluctuations strongly enhanced in reduced dimensionality. Motivated by this unique scenario 
and by the very recent experimental realization of trapped quasi-two-dimensional fermionic clouds, 
we include one-loop Gaussian fluctuations in the theoretical description of resonant Fermi superflu¬ 
ids in two dimensions demonstrating that first sound, second sound and also critical temperature 
are strongly renormalized, away from their mean-field values. In particular, we prove that in the 
intermediate and strong coupling regimes these quantities are radically different when Gaussian 
fluctuations are taken into account. Our one-loop theory shows good agreement with very recent 
experimental data on the Berezinskii-Kosterlitz-Thouless critical temperature [Phys. Rev. Lett. 
115 , 010401 (2015)] and on the first sound velocity, giving novel predictions for the second sound 
as a function of interaction strength and temperature, open for experimental verification. 

PACS numbers: 03.75.Ss 05.70.Fh 03.70.-|-k 


Introduction .—Quantum fluctuations play a crucial role 
in low-dimensional systems, rendering the finite temper¬ 
ature properties of a two-dimensional Fermi gas across 
the BCS-BEC crossover substantially different from its 
three-dimensional counterpart. In particular, in accor¬ 
dance with the Mermin-Wagner-Hohenberg theoreniir— 
for d < 2 there can not be a finite condensate density 
at finite temperature, as the fluctuations destroy the off- 
diagonal long-range order; nonetheless two-dimensional 
systems can exhibit algebraic off-diagonal long-range or¬ 
der, allowing for the existance of a quasi-condensante up 
to a certain critical temperature, due to the Berezinskii- 
Kosterlitz-Thouless (BKT) mechanism^. 

Along with the appearance of algebraic long-range or¬ 
der, as observed for the first time in superfiuid “^He, then 
in an ultracold Bose gas^, in an exciton-polariton gasi 
and very recently in an ultracold Fermi gas^, the other 
fundamental signature of the BKT mechanism at work 
is the universal jump in the superfluid density^, going 
discontinuously from a finite value to zero at the criti¬ 
cal temperature, as observed in thin "‘He film&i^. This 
scenario suggests that in a two-dimensional system the 
role of quantum fluctuations should be crucial in de¬ 
scribing several aspects of the systen>ii, as opposed to 
the 3D case for which one could expect from a mean- 
field theory at least qualitative agreement. Recently the 
strongly-interacting Fermi gas has been the object of nu¬ 
merous Montecarloi^ii^ and experimentali^^— investiga¬ 
tions, and in fact it has been observed that Gaussian 
fluctuations strongly modify both the chemical potential 
and the pairing parameter, particularly in the intermedi¬ 


ate and strong coupling regionsii; it has also been shown 
that the correct composite-boson limit is recovered by 
introducing Gaussian fluctuations^^. 

The determination of a full one-loop Gaussian-level 
equation of state needs, however, a proper regulariza¬ 
tion scheme to remove divergences. In the present Let¬ 
ter we use convergence factors in the pair-fluctuation 
propagato r^b^^ to numerically calculate the T = 0 state 
equation for a system of interacting fermions across the 
BCS-BEC crossover. The aim of the present Letter is the 
investigation of beyond mean-field effects at finite tem¬ 
perature: we calculate the first and second sound veloci¬ 
ties, as a function of the temperature and of the binding 
energy, and then calculate the BKT critical temperature 
from the Kosterlitz-Nelson condition^. In particular the 
predictions regarding the second sound velocity provide 
a benchmark for future experimental investigations: we 
expect it to be open to experimental verification quite 
soon, given the rapid advancements in the realization 
and manipulation of ultracold quasi-2D Eermi gases^S. 
On the other hand the theoretical predictions regarding 
the BKT critical temperature and the first sound ve¬ 
locity are compared with recently obtained experimental 
results^i^, showing good agreement in the intermediate 
and BEG regimes. 

Theoretical framework .—The partition function of a 
system of ultracold, dilute, interacting spin 1/2 fermions 
in 2D, contained in a two-dimensional volume at tem¬ 
perature T, with chemical potential pi can be described 
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within the path-integral formalis m^^i^^ as: 

Z = J ( 1 ) 


the full analytical expression of the inverse pair fluctua¬ 
tion propagator M((5) is reported, along with a more de¬ 
tailed derivation, in the Appendix [BJ The collective exci¬ 
tations at T = 0 are gapless due to Goldstone theorent^S: 


with the following (Euclidean) Lagrangian density 





i’a + gip^ '04- 0t > 


( 2 ) 


where ip a (r, t) and iprr (r, t) are complex Grassmann 
fields, a =t, -j. is the spin index, m is the mass of a 
fermion, having defined /3 = (kBT)~^, ks being the 
Boltzmann constant. The strength of the attractive s- 
wave potential is g < 0, which can be implicitely related 
to the bound state energy2ii2^: 


1 _ 1 ^ 1 
9 2L2 ^ 
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with Ek = fi^k‘^/{2m). In 2D, as opposed to the 3D case, 
a bound state exists even for arbitrarily weak interac¬ 
tions, making eb & good variable to describe the whole 
BGS-BEG crossover. The quartic interaction can be de¬ 
coupled by using a Hubbard-Stratonovich transformation 
in the Cooper channel, introducing in the process the new 
auxiliary pairing fields A(r,T), A(r,They corre¬ 
spond to a Cooper pair, being conjugate to two electron 
creation/annihilation operators^^. Moreover the newly 
introduced pairing fields can be split into a uniform, con¬ 
stant saddle-point value Aq and the fluctuations around 
this value as follows: 


A(r,T) = Ao -f 7 ?(r,T) . (4) 

Neglecting the fluctuation fields rj, fj gives us a sim¬ 
ple mean-field (ME) theory, which is generally unreliable, 
due to the fundamental role of fluctuations in two dimen¬ 
sions, but still constitutes the starting point for more re¬ 
fined approaches. In the T = 0 limit the functional inte¬ 
gral and the k-integrations defining the partition function 
at mean-field level are elementary, so that one finally gets 
the mean-field contribution to the equation of stated!, as 
derived in Appendix 


™ r 2 -r 

n^f{9) = -^{9 + -eBr . (5) 

and the single particle excitations of the mean-field the¬ 
ory are: 


Ek = \J (cfc - + ^0 ■ (6) 

On the other hand one can extend the theory including 
the fluctuations fields at a Gaussian leve l^^i^^ . The re¬ 
sulting equation of state reads: 

, ( 7 ) 
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At finite temperature a gap will appear, however it is 
extremely small below Tbkt as shown in Refi^. Erom 
Eq. (O and Eq. (O we get the full one-loop equation of 
state, by also using Eq. dS]) we get ^ as a function of the 
crossover. 

In conclusion of the present section we note the grand- 
potential in Eq. © cannot be evaluated as is, being 
affected by divergencies related to the modeling of the 
interaction using a contact pseudo-potential rather than 
a realistic one. Many different regularization approaches 
can be used, like the dimensional regularization in 2D 
in the BEG limit^^, the counterterms regularization^ or 
regularization with convergence factor o^^d^ . The first 
two are more suited to obtain analytical results, partic¬ 
ularly in the BEG limit, while the last method has been 
shown to be suited to obtain numerical results across the 
whole crossover—. Wanting to investigate numerically 
the whole crossover, our grand potential is regularized 
by introducing convergence factor s^^d^ : 


r2g(/i) 
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First and second sound .—The first sound velocity 
is calculated from the regularized grand potential in 
Eq. (|n]), by using the zero-temperature thermodynamic 
relation^: 



n / 1 

m\L‘^ d^j? J 
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Using the mean-field equation of state to calculate the 
chemical potential, one would find Cs[iimf) = vfI'/2 
across the whole BCS-BEC crossover, vp being the Fermi 
velocit} d^d^ . Our equation of state with Gaussian fluc¬ 
tuations yields, as expected, a critically different Cg'. it 
slowly tends to the aforementioned value in the BCS 
limit, showing, on the other hand, a remarkable differ¬ 
ence in the intermediate and BEG regimes, tending to 
the composite boson limit derived in Ref—. We plot 
this result in Fig. [1] noting that it exhibits good agree¬ 
ment with the experimental data in Ref—. By adapting 
the thermodynamic approach of Ref— we verified that 
the T-dependence of Cg in the superfluid phase is very 
weak, see the inset of Fig. [T] 

Beside the first sound, propagating through density 
waves, a superfluid can also sustain the second sound, 
a purely quantum-mechanical phenomenon propagating 
through a temperature wave^^. In order to calculate 
the second sound velocity we follow the treatment in2^ 
starting from the free energy of the system, substantially 
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FIG. 1: The first sound velocity at T = 0, calculated using ^ 
and Ao from the Gaussian-level equation of state (black solid 
line), and using their mean-field counterparts (blue dashed 
line), which give a constant value Cs/vf = In the 

strong coupling regime a full Gaussian-level equation of state 
is needed to correctly describe thermodynamic quantities, 
there our prediction correctly tends to the composite boson 
limit (red dotted line). Blue circles and red squares are exper¬ 
imental observations of the sound velocity, as reported in^, 
respectively obtained by directly measuring the speed of a 
density wave and through the equation of state. Inset: tem¬ 
perature dependence for log(es/eF) = —10, —5, 0, 5,10 (from 
top to bottom). 

treating it as a gas of independent single particle and col¬ 
lective excitations, neglecting hybridization through Lan¬ 
dau damping; this approach will be justified shortly when 
discussing the BKT critical temperature. We find the 
fermion single particle contribution to the free energy^: 

and the bosonic one, from collective excitations: 

= ( 12 ) 

^ q 

The total free energy is then F = F^+Fcoi+Fsp where the 
zero-temperature energy Fq is a T-independent constant, 
unimportant as far as the present Letter is concerned. 
The entropy is readily calculated as S' = —{dF/dT)jq^L 2 
and introducing the entropy per particle S = S/N the 
second sound velocity is^Zr— : 



where rig and n„ are the superfluid density and the nor¬ 
mal fluid density, respectively. In contrast with the 3D 
caseZZ here the second sound has a discontinuity at the 
critical temperature, as a consequence of the universal 



FIG. 2: The second sound velocity, as a function of the tem¬ 
perature T/Tf, for varying values of ^bI^f- The character¬ 
istic structure with a minimum followed by a linear increase 
evolves into a constant second sound velocity approaching the 
BEG regime. 

jump in the superfluid density, as also noted inZZiZi; the 
critical temperature will be calculated in the next section. 
Our results are reported in Fig. [21 we note that the sec¬ 
ond sound velocity shows a characteristic minimum in the 
BCS and intermediate regimes, as also noted in the 3D 
unitary caseZZ^ evolving into an approximately constant 
second sound velocity approaching the BEG regime. 

Critical temperature: the Berezinskii-Kosterlitz- 

Thouless transition .—As mentioned in the Introduction 
the low-temperature physics of a 2D attractive Fermi 
gas is essentially different from that of a 3D gas: the 
Mermin-Wagner-Hohenberg theoremZ^— prohibits the 
symmetry breaking at finite temperatures, so that 
one can find off-diagonal long-range order and a finite 
condensate density only at T = Tc = 0. However 
quasi-condensation, i.e. the the algebraic decay of the 
phase correlator (exp(i0(r)) exp(i0(O))) ~ |r|~^ where 
77 is a T-dependent exponent and 9 is the phase of the 
order parameter, is observed up to a finite temperature 
Tbkt, known as the Berezinskii-Kosterlitz-Thouless 
(BKT) critical temperatureZiZ. The other fundamental 
signature of the BKT mechanism is the universal jump 
in superfluid density at the critical temperature, i.e. 

^ ~ 0. The transition temperature is 

determined through the Kosterlitz-NelsonZ condition 

ksTsKT = -^^naiTsKT) (14) 

which allows one to calculate Tbkt, known the super¬ 
fluid density. Within the present framework we write 
the superfluid density as Ug = n — Unj — Un.b where n 
is the density of the system and Unj and Un.b are nor¬ 
mal density contributions arising, respectively, from the 
single particle excitations and from the bosonic collec- 
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tive excitations. Using Landau’s quasiparticle excitations 
formula^ for fermionic: 


J (eoa + 1)2 (15) 


and for bosonic excitations: 


'^n,b — 2 


P f 




(27r)2^ [eP^q-lf 


(16) 


T he single particle e xcitation spectrum is Ek = 
\J{k‘^1 {2m) — /i)^ + Aq, as derived in Eq. ©, the col¬ 
lective excitations spectrum in Eq. dS]) can replaced by 
a phonon-like linear mode ojq ~ Csq within a very good 
approximation as far as the critical temperature is con¬ 
cerned, as we verified. 

As noted in Ref.— Eq. ITSl) and Eq. (fT6l) hold as long 
as there is no Landau damping hybridizing the collective 
modes with fermionic single-particle excitations, other¬ 
wise the bosonic normal density would need to be modi¬ 
fied. In our case one can easily verify that for es > 1 the 
condition cb ^ ksT holds in the whole temperature re¬ 
gion of interest, strongly suppressing the pair breakup 
and the Landau damping^. On the other hand, for 
^ 1 we verify that the critical temperature is deter¬ 
mined by the fermionic contribution to the normal den¬ 
sity, as one would expect, making eventual corrections to 
Un^b neglectable. We then conclude that Eq. dig and 
Eq. (fTHl) correctly describe the normal density for the 
entire superfluid phase. 

By numerically solving Eq. m we find the transition 
temperature Tbkt at different points of the BCS-BEC 
crossover. Our results, shown in Eig. [31 are compared 
with very recently obtained experimental data in Refi^, 
showing a excellent agreement with experimental data at 
least for cbI^f ^ 1 - 

We stress that with respect to other derivations of 
Tbkt in the 2D BCS-BEC crossove r^^d"*^~ — the present 
theoretical prediction of Tbkt includes the contribu¬ 
tion from a Gaussian-level equation of state along with 
the contribution from the bosonic collective excitations. 
These contributions are critical in correctly fitting ex¬ 
perimental data, as clear from in Eig. 0 We find that a 
theory of fermionic only excitations, like those developed 
i n^9i44 Qj. in a slightly different context in^, overestimates 
the critical temperature in the intermediate and strong¬ 
coupling regimes. Conversely, not using a Gaussian equa¬ 
tion of state underestimates the critical temperature in 
the BCS regime, see Eig. [3] 

Moving towards the BGS side of the crossover, how¬ 
ever, the agreement is slightly worse, the experimen¬ 
tal Tbkt being bigger than 0.125eF; by inserting into 
Eq. (fTTl) the relation n = m/{fi^7T)eF it is easily seen 
that the critical temperature is not allowed to exceed 
the value Tbkt = 0.1256^. Thus we conclude that the 
slightly worse compatibility observed cannot be repro¬ 
duced within the framework of the Kosterlitz-Nelson cri¬ 
terion, as defined by Eq. (O, and should be attributed 



FIG. 3: Our theoretical prediction for Tbkt (black solid line) 
as compared to recent experimental observation reported in^, 
temperature estimated through algebraic decay, the error bars 
account for statistical and systematic errors. Our prediction 
uses a Gaussian equation of state, including the contribution 
from single particle modes and collective excitations. A the¬ 
ory with fermionic only excitations (gray dot-dashed line) fails 
to provide an agreement with experimental data in the BEG 
regime, whereas a theory using a mean-field equation of state 
would underestimate Tbkt in the weak coupling regime (gray 
dashed line). The strong coupling extrapolation (red dotted 
line) in Eq. (|17l) gives good results in the strong and inter¬ 
mediate coupling regimes. 


to different physics, like the mesoscopic effects mentioned 
iii2S in the same regime. Nonetheless we stress that our 
results are still within 1 . 2 cr from experimental data, when 
statistical and systematic errors are taken into account. 

Strong coupling limit —In the strong coupling limit an 
attractive Fermi gas maps into a Bose gas, in particular 
by using the relation between the fermionic and bosonic 
scattering lengths ob = ap /and the relation 
cb = Ah ?one finds cbIsf = ^/{nBa'^) with 
K = exp (—27 — I/ 2 )/ 7 r fv 0.061 and ub = np/^. The 
strong coupling regime for the Fermi gas corresponds to 
the extremely dilute limit for the Bose gas. In this limit 
the fermionic contribution to the normal density in Eq. 
(ITSll is neglectable because the energy gap Aq becomes 
extremely large. Moreover the integration for in Eq. 
(fT6)l is analytic, we solve the Kosterlitz-Nelson condition 
and expand in powers of Cs, obtaining the following ana¬ 
lytical estimate for the critical temperature in the strong¬ 
coupling regime: 


kpTpKT 


Tb^f _ 234 

</WS) 3 ( 3C ( 3 ))§ 


(17) 


with ^(3) ~ 1.202 and fiB = f-F/ {nBO?B)) '^4 
Ref—, we report this result in Fig. [31 An alternative 
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estimate of Tbkt can be derived mapping the attractive 
Fermi gas to a Bose gas, we report the derivation in Ap¬ 
pendix [cl noting however that his validity is limited to 
the strong coupling regime. 

Conclusions .—There are several open problems for the 
physics of ultracold atoms which can be faced employ¬ 
ing one-loop Gaussian fluctuations. Here we have shown 
that the Berezinsky-Kosterlitz-Thouless critical temper¬ 
ature of the superfluid-normal phase transition can be 
extracted from an description of the superfluid density, 
which takes into account Gaussian fluctuations in the 
finite-temperature equation of state. The agreement 
with very recent experimental data for both the criti¬ 
cal temperature^ and the sound velocity^! is remarkably 
good and crucially depends on the inclusion of quan¬ 
tum and thermal Gaussian fluctuations. More gener¬ 
ally, Gaussian contributions to the equation of state are 
relevant for Bose-Fermi mixtures^, for unbalanced su¬ 
perfluid fermions^, and also to investigate the conden¬ 
sate fraction in the BGS-BEG crossover—. Finally, we 
stress that in addition to ultracold atomic gases there 
are several other superfluid quantum many-body systems 
where the methods of functional integration and Gaus¬ 
sian fluctuations play a relevant role to achieve a mean¬ 
ingful and reliable theoretical description. Among them 
we quote high-Tc superconductors^, neutron matter in 
the BGS-BEG crossover—, quark-gluon plasma—, quark 
matter in stars^, and, more generally, quantum fluids of 
light^. In particular, our results pave the way for a bet¬ 
ter understanding of the strong-coupling limit of other 
two-dimensional systems with BGS pairing, e.g. bilay¬ 
ers of fermionic polar molecules^ 2 iS or exciton-polariton 
condensates^. 
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Appendix A: Mean-fleld treatment 


Starting from the same Lagrangian as in the main text: 




M 

2m 


■0^ + 5 V't (Al) 


with the same notation, after the Hubbard-Stratonovich 
transformation one obtains the new (Euclidean) La¬ 
grangian density: 


= iio 


fidr - M 

2 m 


i/’t i/’t+Ai/’t i/it- 


9 

(A2) 


and the functional integration needs to be extended over 
A, A. As mentioned in the main text the newly intro¬ 
duced pairing field can be split into a uniform, constant 
saddle-point value Aq and the fluctuations around this 
value as follows: 


A(r,T) = Ao-f r 7 (r,T) , (A3) 

The mean-field approximation consists in neglecting the 
the fluctuation fields 77 , in this case the functional in¬ 
tegral defining the partition function can be carried out 
exactly, yielding: 


= exp |-^^| = exp{-/3Hm/} , (A4) 

where 

%^ = -Tr[ln(Go-i)]-/3L^^ = 
h g 

= - ^ [21n(2cosh(^A«p(A:)/2)) - /S^k] - (A5) 

k ^ 

with — the trace being taken in reciprocal space 

and in the Nambu-Gor’kov space, with 


- 1-1 


G- = 


fi-dr + + /i 


(A 6 ) 


and the single-particle excitation spectrum is found solv¬ 
ing for the poles of the Nambu-Gor’kov Green’s function 
Go in momentum space^^: 


Ek — \J (cfc ~ + Ag . (A7) 


In the T = 0 limit the k-integration in Eq. (IA5I1 can be 
carried out analytically in the 2D case, one gets: 


Aq) — 


mL^ 
47rft^ L' 


9 ^ + MV + Aq + -Ag -I- 


-Aoln ( 


—g + + Ag 


£b 


.(A 8 ) 


Imposing the saddle-point condition for Aq, i.e. 
{dflmf/dAQ)^ y = 0 , one obtains the gap equation: 


Aq — \l 2eb{g + -^^b) 


(A9) 


plugging this result back into the MF grand potential we 
get the MF equation of state: 


^ . 1 ,9 

f^m/(M) = -^(m+2^3) ■ 


(AlO) 
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Appendix B: Gaussian fluctuations 

Restoring the fluctuation fields 77 , ^ at a Gaussian level, 
the partition function readsi^: 


As noted in the main text in this limit an attractive Fermi 
gas maps to a Bose gas, each boson has mass ms = 2m f 
and the density is u_b = uf/S. The Fermi energy of the 
original Fermi gas is then: 


Z = J VrjVfj exp , (Bl) 


h^TT h^TT 

ep = - np = 4 - ns 

mp m_B 


(Cl) 


where 

= \Y.^n{Q),ri{-Q)) M(Q) ) (B2) 

Q 

having introduced the Fourier-transformed version of the 
fluctuation fields, with Q = (in„,q), = 2nnl(3 being 

the Bose Matsubara frequencies. The 2x2 matrix in Eq. 
(IB2II is the inverse propagator for the pair fluctuations, 
its matrix elements are defined hy^: 


From now to the end of the present section F and B sub¬ 
scripts will be used to distinguish between, respectively, 
fermionic and bosonic masses, densities and scattering 
lengths. The binding energy ep and the fermionic scat¬ 
tering length ap are related by the equation found by 
Mora and Castin^: 


e^'^ma?p 


(C2) 
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= E' 


E 

k 


-E-E' 

1 


E- 

1 


E' 


+ E + E' iuJn — E — E' 


(B3) 

|(B4) 


where u = uw = ^u = Uk = V'l-^k, 

u — lAk-fqi ^ B — F/sp(k), E — F/sp(k-t-q). The 

remaining matrix elements are defined by the relations: 
^ 22 ( 9 ) = Mii(- 9 ), M[ 2 i(g) = Mi 2 (g). The quantity l/g 
appearing in the definition of Mu (Q) is removed by using 
the scattering theory result: 


1 1 1 

9 2Z/^ k 


(B5) 


By integrating out the 77 (r,T), 7 ?(r,T) fields in Eq. (IBII) 
we get the Gaussian contribution to the grand potential: 

Ao) = ;^^ln(detM(g)) (B6) 

g 

and the equation of state {n) is found, like in the ME 
case, by inserting the gap equation from Eq. dMl). By 
imposing the condition det(M) = 0 one can find the col¬ 
lective excitation spectrum, which will have, in the low- 
momentum limit, the following expression 


It can combined with the relation found between 
the bosonic and fermionic scattering lengths^ ap = 
apj to obtain: 


££ — ^ 
ep npa?p 


(C3) 


where k = exp (—27 — l/ 2 )/ 7 r Ri 0.061, 7 being the 
Euler-Mascheroni constant. The Berezinskii-Kosterlitz- 
Thouless critical temperature for a dilute 2D Bose gas^ 
has been estimated using Montecarlo techniques^SiSi: 


TpKT _27r_ 

np mp\og{^^^) 


where 


UeS 


Att 

mp log(l/na|) 


(C4) 


(G5) 


and Montecarlo simulations yield ^ ^ 380. Moreover Eq. 
(IC4I) can be rewritten by using Eq. (IGip and Eq. (IC3I) 
as: 


TpKT _ 1_1_ 

2log(ilog(K-if|)) 


(C 6 ) 


hUq = €q (Acg -I- 2mcl) (B7) 

A and Cg being a function of the crossover. As already 
noted in the main text the collective excitation spectrum 
is gapless as a consequence of Goldstone theorem. 

Appendix C: Alternative determination of the 
critical temperature in the Bose limit 

An alternative estimate of the BKT transition temper¬ 
ature can be given by extending to the intermediate cou¬ 
pling regime the strong coupling result valid for ^ ep. 


This result can be compared with the experimental data 
reported in^; we observe that the theoretical prediction 
in Eq. (IG 6 I) correctly fits experimental data within the 
reported statistical errors. Gonversely we can leave ^ as a 
free parameter and estimate it through Eq. (IC 6 I) and the 
experimental data using a simple least squares method, 
the result is ^ = 554 ±179 which is compatible with 
the Montecarlo estimate in^. However we must stress 
that this alternative result in Eq. (IC 6 I) is essentially a 
composite-boson extrapolation from the strong coupling 
regime which is not reliable as we approach the BGS side 
of the resonance, as demonstrated by the divergence in 
TpKT when epjep = exp( 47 r/^)«: ^ 0.063. 
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FIG. 4: The condensate fraction in the BCS and interme¬ 
diate regimes, as calculated from Eq. (ID2II . compared with 
Montecarlo data from Ref.—. 

Appendix D: Condensate fraction 

A fundamental quantity in the study of ultracold sys¬ 
tems is the condensate fraction. In a 2D system one ex¬ 
pects a finite condensate density Uq only a T = 0 due to 
the Mermin-Wagner-Hohenbergii^ theorem. In the case 
of a 2D attractive Fermi gas the condensate density is 
given by^ 

no = lim G2i{p,ium)Gi2{p,iuin) (Dl) 

/3^oo ^ ' 

p,m,n 


where G is the single-particle Green’s function, oj„, ujm 
are Fermi Matsubara frequencies, and the condensate 
fraction is simply nQjn. In the BCS limit all the integrals 
can be carried out analytically as in Ref—, yielding 



We use the Gaussian equation of state for p and Aq to 
compare this theoretical prediction for the condensate 
fraction with the recent Montecarlo results in Ref—, as 
reported in Fig. 01 in the BCS and intermediate regimes. 
Clearly the rather good agreement breaks as the inter¬ 
action gets stronger; an extension of this analysis to the 
BEC regime would require the evaluation of Eq. m 
across the crossover, recalling that for self-consistency 
the single-particle Green’s function is to be calculated at 
one-loop level^. This calculation is the subject of ongo¬ 
ing work by the authors. 
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